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Abstract 

1 

The  effective  velocity  of  esoape  of  the  powderl  gas  from  a  gun,  as 
defined  by  Bq»  (l) ,  Is  a  quantity  that  is  independent  of  the  mass  of  the 
recoiling  parts*  A  knowledge  of  this  esoape  velooity  permits  the  calcula¬ 
tion  of  the  maximum  velooity  of  free  recoil,  a  quantity  of  interest  in  the 
theory  of  damped  recoil,  and  therefore  of  value  in  the  design  of  a  reooil 
mechanism*  Although  this  esoape  velooity  is  independent  of  the  reooiling 
mass,  it  depends  on  other  faoters,  so  that  the  oomujai  assumption  vis*,  that 
it  has  a  universal  value,  is  unsatisfactory* 

It  is  the  purpose  of  this  report  to  develop  a  j  method  for  predicting 
the  esoape  velooity  that  will  be  satisfactory  for  values,  of  the  ratio  of 
the  mass  of  the  charge  to  that  of  the  projectile  as  large  as  unity*  The 
theory  follows  Hugoniot  in  treating  the  efflux  of  gae  After  departure  of 
the  projectile  as  equivalent  to  the  emptying  of  a  reservoir  through  a  nozzle, 
values  of  quantities  at  the  breech  being  used  for  values  in  the  reservoir* 

In  order,  however,  that  the  theory  may  hold  for  high  values  of  s  ,  values 
at  ejection  of  the  ratios  of  mean  pressure  and  mean  density  to  breeoh  values 
are  derived  from  the  Pidduoto»Ksnt  special  solution  for  the  motion  of  the 
powder  g&s* 


The  resulting  equation  is  applied  to  the  240mm  howitzer.  Model  1918,  to 
the  3”  eeacoast  gun.  Model  1898,  and  to  small  arms*  The  agreement  for  the 
large  caliber  guns  is  within  3#  on  the  average,  the  value  of  e  being  as  large 
as  l/5  for  the  3"  gun*  It  is  shown,  moreover,  that  the  extension  to  valueB 
of  s  as  large  as  unity  does  not  involve  more  than  about  a  5$  correction;  thus 
even  if  the  latter  should  be  in  error  by  20$,  only  1$  possible  error  would  be 
added*  The  equation  is  thus  expected  to  be  equally  satisfactory  for  6  =  1. 

A  complete  list  of  the  main  equation  and  of  the  auxiliary  formulas  is  given 
in  the  summary* 


As  by-products  there  are  developed  formulas  for  the  explicit  calculation 
of  the  main  parameter  &_  of  the  Kent  solution  as  a  function  of  e  and  a 

1+l4for 

+  e/J 


nunmrioal  verification  that  the  value  valid  for  small  &  ,  viz* 


the  ratio  of  mean  pressure  to  breeoh  pressure  before  ejection  is  still  ao curate 
even  for  rather  large  values  of  e  •  These  results  should  be  useful  in  ordinary 
pre-ejection  interior  belli at ios* 
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The  design  of  a  system  for  damping  the  recoil  of  a 
gun  requires  a  Knowledge  of  the  maximum  velocity  that  the  gun 
would  acquire  if  allowed  to  recoil  with  no  retarding  force 
whatever.  The  latter  quantity,  called  the  maximum  velocity  of 
free  recoil,  we  shill  denote  by  the  symbol  Vr^»  (The  subscript 

X 

r  denotes  "recoil'1  and  the  subscript  f  "final  value".)  It  is 
commonly  determined  by  one  of  two  methods.  The  first  consists 
in  mounting  the  gun  on  rollers,  designed  to  be  as  friction-free 
as  possible,  and  determining  directly  the  actual  velocity  of 
recoil  as  a  function  of  time.  For  this  purpose  some  form  of 
velocimetbr  or  chronograph  is  used.  The  second  method  consists 
in  mounting  the  gun  as  a  ballistic  pendulum  and  determining  the 
rise  in  the  center  of  gravity  of  the  whole  System  when  the  gun 
is  fired.  From  this  measurement  the  maximum  velocity  is  easily 
calculated  for  the  free  recoil  of  the  compound  system.  To  deter¬ 
mine  the  maximum  velocity  of  free  recoil  of  the  gun  alone  one 
then  assumes  that  the  recoil  momentum  is  independent  of  the  mass 
of  the  recoiling  parts.  Such  an  assumption  is  equivalent  to  the 
assumption  that  the  gas  pressure  inside  the  gun,  as  a  function 
of  space  and  time,  is  independent  of  the  mass  of  the  recoiling 
parts.  The  latter  assumption  is  entirely  acceptable,  in  view  of 
the  known  smallness  of  the  effect  of  recoil  on  interior  ballistic 
processes.  If  Mp  denotes  the  total  mass  that  recoils  in  the 

pendulum  experiment  and  Mr  denotes  the  mass  of  the  recoiling  parts 

of  the  gun  in  an  ordinary  firing,  the  desired  value  of  Is 

“r 

therefore  given  by  the  product  of  the  pendulum  velocity  and  the 
factor  Mg/Mv.  The  roller  method  is  used  only  for  large  caliber 

guns,  while  the  pendulum  method  is  practicable  only  for  small 
calibers,  although  it  has  been  used  for  guns  as  large  as  the 
three-inch. 


The  results  of  such  experiments  are  commonly  presented 
by  giving  the  value  of  a  quantity  called  "the  effective  velocity 
of  escape  of  the  powder  gas",  denoted  by  V  .  This  quantity  is 

defined  by  the  following  equation: 

Vi*  ’  Mvm  +  cve*  « 

± 

where : 


-3- 


M  =  mass  of  recoiling  parts 
K  £.  mass  of  projectile 


C  =  mass  of  powder  charge 


V"  *.=  maximum  velocity  of  free  recoil 

n 

=  muzzle  velocity  of  projectile 


£■  inc  e  Eq . 
we  may  interpre 


(1)  obviously  refers  to  a  conservation  of  momentum, 

t  V  as  follows,  ihe  final  momentum  M  \  of  the 
e  1  it 


freely  recoiling  gun  must  be  equal  to  the  sum  of  the  momentum  MV 
of  the  projectile,  the  final  momentum  of  the  escaping  powder  131 
gas,  and  the  momentum  imparted  to  the  surrounding  air  by  the 
blast.  Thus  CV  represents  the  sum  of  the  last  two  momenta.  We 
now  make  the  reasonable  assumption  that  A  and  MV  are  independent 


oi'  the  state  of  the  atmosphere,  so  that  their  values  would  be 
unchanged  if  the  firing  were  into  a  vacuum.  It  follows  that  V 
represents  the  final  value  that  would  be  achieved  on  the  average 
by  the  axial  component  of  powder  gas  velocity,  if  the  firing 
were  into  a  vacuum.  The  quantity  V.  is  not  the  same  as  the 

average  gas  velocity  at  the  muzzle  after  ejection  of  the  projectile, 
since  the  powder  gas  speeds  up  on  leaving  the  gun,  just  as  would 
the  gas  in  a  nozzle  in  passing  from  the  convergent  to  the  divergent 

portion.  Furthermore  we  shall  calculate  7^  not  by  ary  direct 
consideration  of  processes  at  the  muzzle,  '  “ 

of  V  ,  arid  use  of  Eq.  (1).  For  these  reasons,  although  we  have 
interpreted  Eq.  (1)  as  referring  to  a  conservation  of  momentum, 

it  appears  desirable  to  adhere  to  the  position  taken  above,  viz. 
that  Eq.  (l)  is  simply  a  definition  of  V^. 

The  importance  of  the  quantity  V  lies  in  the  fact  that  it 

is  presumably  independent  of  the  mass  Mr  of  the  recoiling  parts. 

One  can  understand  this  statement  in  the  light  of  the  following 
remarks.  In  the  usual  case  where  M  >:>  M,  it  is  known  that  recoil 
has  little  effect  on  processes  insiSe  the  gun.  Thus  gas  pressure, 
as  a  function  of  space  and  time,  will  depend  only  very  weakly 

on  the  ratio  M/M  .  so  that  the  momenta  MrV  ~  and  MV  will  not 

r 1  *  ix  m 

depend  appreciably  on  Mr.  Then  from  ;-q.  (1)  the  escape  velocity 
V  will  be  independent  of  m  .  The  lack  of  influence  of  M  on 

S  f  J- 

Vg  makes  the  latter  a  useful  quantity  for  the  expression  of  results 

on  recoil,  since  various  values  of  A  may  be  contemplated  in  a 

design.  Although  V  is  thus  independent  of  M  ,  we  shall  see, 

however,  that  it  may  depend  on  other  variables,  so  that  the 
common  practice  in  ordnance  engineering  of  using  a  universal 


-ht- 


value*  for  V  is  therefore  not  justified. 

3 •  The  Quantity  K 
We  define  K  by: 

K  s  Ve/Vm  (2) 

Then,  from  (l)  and  (2) : 

Mrvr$  ■  <M  +  KC)\  (3) 

II.  Theory  of  the  Velocity  of  Free  Recoil 

1.  Ihe  Velocity  of  Free  Recoil  Before  Ejection 

Let  Vrni  be  the  velocity  of  the  recoiling  parts  when  the 

base  of  the  projectile  is  flush  with  the  muzzle.  For  values  of 
e  =  C/M  that  are  not  too  large  (e.g.  e<l/3),  all  theories  of  the 
motion  of  the  powder  gas  give  the  value  1/2  CVm  as  the  momentum 

of  the  powder  gas.  The  recoil  momentum  at  this  stage  is  then  given 
by: 

Vn  *  (M  +  2  c)Vfflj  U) 


2 .  Additional  Velocity  of  Free  Recoil  After  Ejection 

A.  Discussion  of  the  Problem 

After  the  projectile  has  left  the  gun,  there  are 
still  several  forces  acting  that  tend  to  change  the  momentum  of 
the  freely  recoiling  gun.  These  are 

(a)  The  force  due  to  gas  pressure  on  the  breech. 

(b)  The  force  due  to  gas  pressure  on  the  muzzle  ring, 

i.e.  the  force  due  to  the  pressure  of  escaping 
gas  acting  on  an  area  J(L|-D£),  where  D2  is  the 

, outride  diameter  at  the  muzzle  and  D-,  the 
inside  diameter.  1 

(c)  The  drag  due  to  gas  friction  on  the  walls  of  the 

chamber  end  the  bore. 

(d)  The  chambrage  force,  i.e.*  the  force  due  to  gas 

pressure  acting  on  the  cpne  joining  the  chamber 
or  cartridge  case  with  the  bore.  Let  Fig.  1 
schematize  roughly  the  chamber  (or  cartridge 
case)  and  the  bore  of  a  gun. 


*  E.g.  the  well  known  texts  by  Tschappat,  MacFarland,  and  Hayes-  all  suggest 
the  value  Vfi  =  4700  ft/sec.  1  B 


nozzle.  In  the  usual  resevoir  problem  one  takes  the  gas  velocity 
to  be  zero  in  the  resevoit;  we  shall  therefore  interpret  resevoir  , 
values  as  referring  to  ,the  breech.  Let  p  an&  p  denote  ttie  resevoir 
or  breech  values  of  pressure  and  density  at  time  t  and  p  and  p 
their  space-mean  values.  The  subscript  e  wij.1  denote  initial  or 
shot  ejection  values.  Let  us  now  define  r^  and  r2  by: 


P  =  rxP 

(6) 

P  =  r2? 

(7) 

We  shall  assume  r-^  and  r2  to  remain  constant  during  efflux  and 


therefore  equal  to  their  shot  ejection  values,  which  can  be  cal¬ 
culated  from  the  theory  of  flow  before  ejection.  With  neglect  of 
corrections  due  to  gas  imperfection,  the  rate  of  mass  efflux  is 
given  by: 


<P  *  A  7 


1/2  / _ 2x 

'Y  +1'' 


■Y-+1 

27r=t) 


(pp)1/2> 


(3) 


where  y  is  the  ratio  of  specific  heats  of  the  powder  gas.  Eq.  (8) 
holds  only  when  the  ratio  of  external  pressure  to  breech  pressure 
is  less  than  a  certain  critical  value,  but  by  the  time  the  breech 
pressure  has  fhllen  so  low  that  this  critical  value  is  exceeded, 
the  gun  is  practically  emptied  of  gas,  so  that  no, appreciable 
error  can  occur  on  this  account.  The  pressure  p  and  density  p 
are  connected  by  the  ddiabatic  relation: 

1 

PP~Y  c  PePe  *  (9) 


If  X  denotes  the  ratio  of  the  total  volume  from  breech  to  muzzle 
to  the  area  A  of  the  bore,  the  differential  equation  for  emptying  is 


a  y  ilfi 
A  *  dt 


-  cp 


(10) 


Combination  of  Eqs. (6) ,  (8),  (9),  and  (10)  and  integration  with 
the  initial  condition  p  =  p  gived'T 

p  =  Pe(i  +  1  >  (11) 

JT»-1. 


where 


T  =  (X/r1)TrxK  Y  \pe/ve)1/2, 


a  time  constant  that  governs  the  rate  of  fall  of  pressure  after 
ejection.  Let  T  denote  the  space-mean  absolute  temperature  of 

the  powder  gas  at  shot  ejection,  and  F,  the  gas  constant  per 

unit  mass.  Then,  from  the  equation  of  state: 


Pe^i^e 

(13) 

Also,  from  the  definition  of  p, 

Pe  =  c/ (ax) 

(14) 

Integration  of  (11),  with  use  of  (12),  (13) ,  and  (14) ,  gives  for 
the  additional  recoil  momentum  after  ejection: 


where 


<XJ 

a  r  pdt 

o 

f(r)  =  v 


-1/2  1/2  1/2 
=  *1  r2  i'CY)CCR1Te) 


(15) 


(16) 


The  usual  uniform  density  theory  for  gas  flow  before  ejection 
gives  r, =1  and  r  =  (l+e/2)/(l+s/3) .  These  values  will  be  valid 

■  "S 

for  values  of  e  =  C/M  that  are  not  too  large,  say  for  e  <  V,re 

shall  consider  later  what  modifications  of  these  values  may  be 

necessary  for  larger  values  of  e. 

C.  Estimation  of  the  Election  lamp era lure 

During  the  travel  of  the  projectile  in  the  gun 
energy  is  lost  from  the  system  (powder  gas  plus  projectile) 
through  two  effects,  viz,  by  heat  transfer  from  the  hot  powder 
gas  to  the  cooler  walls  of  the  gun  and  by  bore  friction  of  the 
projectile.  Let  k  denote  the  ratio  of  the  total  energy. lost  in 
these  ways  (at  the  moment  when  the  projectile  base  has  just 
reached  the  muzzle)  to  the  kinetic  energy  1/2  MV^  of  the  projectile. 
The  kinetic  energy  of  the  powder  gas  is  given  by  l/2  C/b  V^,  where 

h  has  the  value  3  for  sufficiently  small  values  of  e;  its  value 
will  be  given  later  for  larger  values.  The  energy  equation  of 
interior  ballistics  then  gives 

Bpe  =  X  -  nl(l±!±eZ&)  ,  (17) 


where  X  denotes  the  specific  force  of  the  powder.  The  best 
estimate  available  for  k  is  given  by  a  statistical  analysis  ot 
certain  firing  records  for  single-perforated  powder.  Such  an 


analysis  showed  that  k  depends  significantly  on  w,  D>  A,  and  T 
only,  where  ° 


w  =  web  thickness  of  the  powder 
D  =  effective  caliber  of  the  gun  =?\J 
A  m  density  of  loading 

Tq  =  adiabatic  flame  temperature  of  the  powder 
A  least  square  fit  of  log  k  as  a  linear  function  of  the  logarithms 


of  w,  D,  A,  and 


1000 


-  1  gave  coefficients  close  to  +1,  -1,  -1, 


and  -~2  respectively.  The  final  formula  obtained  was: 


(18) 


11000 


In  Eq.  (18)  w  and  D  are  to  be  expressed  in  the  same  units,  A  in 
gm/cm3,  and  Tq  in  degrees  Kelvin,  while  k  is  a  pure  number. 

Eqs.  (17)  and  (18)  then  permit  the  calculation  of  Rjjr  . 

3.  The  Escape  Velocity  for  Small  Values  of  the  Ratio  of 
Charge  to  Projectile  Mass 


Combination  of  Eqs. 


of  the  values  r1  =  1  and  r2 


:  i ; 


(4),  (5),  end  (15) ,  with  use 
)/(l  +  e/3) ,  gives: 


+  e/2 


K  = 


V  /v  =  ^  +  (-1 

e'  m  2  v 


±_eZ2} 
1  ^rT73) 


1/2  1/2 

r(v)(R1Te)  /Vm  , 


(19) 


where  f(-y)  is  given  by  (16)  and  R1T0  by  (17)  and  (18). 

4.  Extension  to  Values  of  e  as  Large  as  Unity 

We  how  consider  what  changes  have  to  be  made  in  (19) 
for  values  of  e  that  are  not  small  compared  ito  1.  For  this 
purpose  wejj-use  thePidduck*-Kent**  solution  |for  the  motion  of 
the  powder  gas  before  ejection.  Reference  is  made  especially 
to  the  treatment  as  developed  by  Kent.  For  convenience  let 


*  Love  and  Pldduck,  Phil.  Trans.  Roy.  Soc.  222.  167(1922) 
**  R.  H.  Kent,  Physics,  2,  319  (1936). 


(20) 


r;'-  -'-a  ■  {.■■■  >) 

Then  the  momentum  of  the  powder  gas  just  before  ejection  is 
given  by: 

Momentum  =  ^  C  Vmh,  (21) 

where  h  is  a  correction  factor  close  to  unity  even  for  e  =  1. 

It  has  the  value: 

-n-1  "I 

l-aQ)  -  lj  ,  (22) 

where  a  is  the  solution  of  the  equation 
o 

2a0(n+l)S(l-a0)_n'1  =  e  ,  (23) 

where  1 

S£/d 

o 

a  function  of  aQ  and  n. 

Expansion  gives,: 

h  =  1  “  §  a0  +  fo(n“6)ao  +  **•  (25) 

We  are  concerned  with  values  of  n  abput  equal  to  4  and  with 
values  of  aQ  less  than  0.1  (for  e  -  1,  aQ  ^  0.1). 

Thus,  to  1  part  in  1090&  ■)  1 

h  =  1  -  §  aQ  (25.1) 

(The  calculation  of  aQ  will  be  considered  ^hortly.) 

Eq.  (4)  then  becomes 

“  (M  +  |  c>  vm  <26> 

It  leads  to  replacement  in  (19)  of  the  first  term  ^  by  h/2. 

We  now  have  to  consider  r,  and  in  (15).  From  Kent's 
solution,  we  derive:  * 


-10- 


(27) 


2n+3 


:(i  + 


With  the  use  of  (23) ,  (24) ,  and  (27)  we  coastruet  the  following 
table  for  the  case  n  =  4  (for  integral  values  of  n,  S  is  a 


polynomial  in  aQ): 


Table  I 


1W3  ; 

1+6/2 

(Value I  of 
l/r2  for  small  a) 


l-e/6 

(Value  of  l/r2 
for  very  small 


0.08 

0.10 

0.20 


0.900725 

0.878106 

0.776940 


1.0933 

1.4871 

4.7421 


0.8788 

0.8520 

0.7361 


0.8822 

0.8578 

0.7655 


0.8178 

0.7522 

0.2906 


As  e  diminishes,  l/r2  approaches  the  value  which  is  labeled 

in  Table  I  as  the  value  for  ’’small  e".  Of  course,  if  e  is 


•siifrie'iMtHS  small  can  be  expressed  ajs  1  -  e/6,  which  is 


accordingly  labeled  in  Table  I  as  the  value  of  l/r2  for  "very  small  e" 

Inspection  of  Table  I  shows  that  the  value  for  ’’very  small  e" 
approximates  l/r2  very  poorly  for  values  of  e  greater  than  unity, 
but  that  the  value  for  ’’small  s”  is  only  4$  in  error  for  a  value 
of  e  as  large  as  4*7  and  only  0.4$  in  error  for  e  -  1.1.  Thus 
for  e  -  1  we  have,  correct  to  better  than  0.2$: 


For  r^  Kent’s  solution  gives: 


r,  *  s  =/(l-V2)n  d|i 


(29) 


Expansion  gives  ^  =  1  -  |  aQ  +  aISt-I).  ae 


af  +  ... 


(30) 


so  that  r’1/2  *  1  -  |  aQ  +  ^(13n  -;18)a£+... 


(31) 


For  n  =  4  arid  e  *  1.1  (so  that  a  530.08),  the  quadratic  term 

amounts  only  to  0.2$.  Thus,  accurately  enough  for  all  values 
of  e  -  1,  we  have; 
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(32) 


>-1/2  - 


-  1  - 


n 


-  ao  =  h,  by  (25.1) 


'1  be  value  of  h  in  (.17)  Is  given  exactly  by 


1 

% 


2n+3 


1  _  2  fn+l') 
e 


a 

*-  o 


We  now  consider  the  calculation  of 
(24)  one  can  obtain  e  as  a  series  in  a 


0 


(33) 


“o'  Frore  Ets-  (23)  arid 
and  can  then  invert  the 


senes 


(by  Taylor's  theorem)  to  obtain  aQ  as  a  series  in  e. 

Unfortunately  this  series  converges  so  slowly  that-  even  if  terms 
are  xept  through  e5,  the  error  amounts  to  for  n  =  A  and  e  =  1  l- 
such  an  error  is  too  great  for  use  of  the  series  in  (33),  V/here  *  * 
accuracy  is  lost  by  subtraction.  Reference  is  made  to  the  Appendix 
for  derivation  of  the  following  procedure.  For  n  =  4  and  e  as 
large  as  1.3  it  gives  correct  to  1  part  in  7000. 


o 


Let 

atn+T) 

=  E 

I 

^ _ / 

ao 

!-ao  ~ 

b 

(35) 

Calculate 

bl  “ 

(1  +  4  nE)1^2  -  :.j 

(36) 

and 

t=J 

III 

r 

15 

(37) 

Th  en 

•>■-!=  1 

(1+  ®  nE1  -  1 

1  (38) 

and 

..  _  b 

Qo  "  l+b 

(39) 

On  replacing  the  1/2  in  (19)  by  h/2.,  as  mentioned  after 
Eq-  (26) ,  and  comparing  (15),  (28),  and  (32)  with  (19),  we  have 
as  our  corrected  formula: 


f('r)(E1ie)1/2/vI„ 


(40) 


where  h  is  given  by  (25-1),  f(_y)  by  (16),  R  T  by  (17)  (18)  , 
and  (33),  and  aQ  by  Eqs.  (34)  to  (39)  .  Any  uncertainty  as  to 

the  validity  of  Eq.  (40)  v/ould  arise  from  the  second,  i-a,  the 
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post- ejection  term.  We  shall  therefore  compare  it  with  ex¬ 
periment  by  multiplying  the  second  term  by  a  "fudge  factor"  F 
and  calculating  in  each  case  the  value  of  F  that  has  to  be 
assigned  to  obtain  agreement  with  experiment*.  We  have  then; 


K 


=  |  +  F  h  ( 


1  + 

1  + 


1/2/V  . 

7  m 


(41) 


'  III.  Comparison  with  Experiment 

1.  Kent’s  yiringg  in  the  2A0mai  Hnwlt^er,  Model  of  1918  Ml. 

These  firings  are  reported  in  Ordnance  Technical  Motes 
No.  6#.  The  following  data  are  relevant;  chamber  volume  =  1790  in3, 

A  =  71.11  in2,  D  =  (A  kp--  9.515  in,  pyro  powder.  Other  data 

7 r 

follow  round  by  round  .1  , , 


r  1 Table  II 


Rd.  N0 • 

w(in) 

e 

*  . 

A,(gm/cm3) 

Vm(f/s) 

Ve(f/s) 

K 

18 

0.036 

0.02825 

0.1546 

822 

4082 

4.966 

19 

0.036 

0.04085 

0.2242 

1021 

3782 

3.704 

9 

0.0438  . 

0.06314 

0.3450 

1056 

4300 

4.072 

11 

0.0488 

0.10453 

0.5751 

1546 

4400 

2.846 

For  y  t  TQ,  and  X  we  use  nominal  values  for  pyro  powder, 

given  by  Hirschf elder,  Kershner,  and  Sherman.**  These  are  Y  =  1.228, 
Tq  =  (2610  +  5/w)  degrees  Kelvin,  and  X=  (303,000+600/w)f t .lbwt/lb, 

where  w  denotes  the  web  thickness  in  inches.  (The  correction  terms  in 
l/w  allow  for  the  effect  of  moisture  and  volatiles) .  We  then  have 
n  =  1/0.228  and  f(v)  s  1.3728.  l'able  III  then  follows,  where  B 


denotes 

the  coefficient  of  F 

in  Eq. 

(41) . 

i 

Table 

Ill 

Rd,  NO. 

k 

ao 

h/2 

1/6 

(R1Te)l/2 

B 

F 

18 

0.4776 

0.00259583 

0.499 

0.3323 

2500  f/s 

4.176 

1.07 

19 

0.3293 

0.00373645 

0.499 

0.3319 

2528 

3.402 

0.94 

9 

0.2932 

0.00573059 

0.498 

0.3312 

2735 

3.559 

1Q0Q+0 . 0. 

11 

0.1759 

O.O0935453 

0.497 

0.3270 

2642 

2.349 

Mean  = 

l.OO+O.o: 

1.00 

^Ordnance  Technical  Notes  No.  6,  "Experiments  in  Interior  Ballistics",  Office 
of  the  Chief  of  Ordnance,  U.S.  Army,  1925. 

**1UDRC  Armor  and  Ordnance  Report  No.  A-204,  page  87. 


The  aalues  given  for  a Q  are  of  course  not  really  accurate  to  so 

many  significant  figures,  but  they  are  accurately  consistent  with 
the  values  of  e  and  v .  Thus  loss  of  accuracy  by  subtraction 
is  avoided  in  Eq.  (33;  for  l/ft.  For  rounds  9  and  11  the  values 
of  V  were  given  only  to  the  nearest  hundred  feet  per  second,  i.e.* 
to  \J2  part  in  44  or  to  1%,  so  that  the  corresponding  values  of  F 
are  uncertain  by  1$.  The  mean  value  of  F  comes  out  exactly  1.00. 


2.  Pastorlza1 s  Firings,  in  the  3”  Seacoast  Gun,  , Model  of  1898 

These  ballistic  pendulum  firings  are  reported  in  the  First 
Progress  Report  on  Ordnance  Board  Program  N6»  2337-2,  "Gas  deflector 
for  Reducing  Recoil  in  Guns",  25,  1^3#*  The  data  follow: 

Projectile  mass  M  *  15.00  lb« 

From  No.  1676,  "Table  of  O.S.  Army  Cannon,  Carriages,  and  Projectiles", 
March  24,  1904,  Revised  Jan*  15,  1924: 

Chamber  volume  *  200  in3 
A  =  7.279  in2 
D  =  3.044  in 


Hirschf elder 1 s 
,4e&rees  Kelvin,  and 
Sec“? .  Various  charges 


No  information  was  given  about  the  powc.er.  All  powders  at 
that  time,  however,  were  pyro  powders  and  the  above  cannon  table 
gives  as  the  typical  web  thickness  w  *  0.047” 
nominal  values  are  then  1=  1.228,  T  =,271 
X  =  316,000  ft .lb.wt/lb.  =  (1017)1  °104fft 
were  fired  without  a  gas  deflector,  but  fiiuzile  velocities  were 
taken  only  on  rounds  16,  17,  and  18  with  a  charge  C  =  5.00  lb. 

The  mean  value  was  Vm  -  1/3(2657  +  2645  +  2660)  -  2654  ft/sec. 

Rounds  14,  15,  19,  and  27  were  the  only  sucdessful*rounds  at  a 
charge  of  5  lb.  without  a  gas  deflector,  bn6  no  muzzle  velocities 
were  taken  on  these  rounds.  We  shall,  however,  use  the  above 
value  V  =  2654  ft/sec  for  these  rounds  also.  Various  recoiling 

masses  were  used.  As  discussed  in  the  introduction,  however,  the 
recoil  momentum  is  expected  to  be  independent  of  Mr,  so  that  it  is 

appropriate  to  use  an  average  value.  The  values  of  ln  sluS 


ft/sec  were  respectively  1970,  1942,  1970,  and  1974,  the  average 
being  1964  slug  ft/sec  or  63,190  Ib/ffc/sec JoB4(l)  then  gives: 


5  Ve  +  15(2654)  =  63,190, 


(42) 


leading  to  V0  =  4676  ft/sec  and  K  =  =  1.762. 

aCaltulatdxaa  then  gives  e  =  1/3,  A  =  0.6921  gm/cm3. 


*  i.e.  successful  in  giving  reliable  values  for  recoil  velocity. 


k  =  0.4393,  aQ  =  0.0276919,  §  =  0.490,  1/h  =  0.3224, 

(r,T  )1^2  =  2538  ft/sec,  B  =  1.318,  and  1  =  0.97.  The  check 
is  again  very  good. 


3.  Small  Arms 


For  small  arms  we  use  some  ballistic  pendulum  firings 
reported  by  Kent  and  Hitchcock.*  These  firings  were  all  with 
IMR  powder  for  which  the  nominal  constants  are  \  =  1.246,  TQ=2795 

degrees  Kelvin,  X=  334,000  ft.lb.wt/lb.,  f  (\)  =  1.3547. 


Da  ta  i 


A.  Caliber  0.50  Machine  Gun  (M19jZl  A. A.') 


w  =  0.022  in 

Vm  =  2500  ft/sec 

V  =  4679  ft/sec 
0 

C  -  245  grains 
M  -  720  grains 
A  =  0.2011  in2 
D  =  0.506  in. 

Chamber  volume  =  1.05  in3 


e  =  0.3403 
A  =  0.922  gm/cm3 
k  =  0.908 
aQ  =  0.0299791 

§1  =  0.490 

1/b  =  0.3224 

(R-,TJ1/2  =  2488  ft/sec 

X  © 

B  =  1.354 


Calculated; 


K  =  1.872 


F,  -  l.o; 

The  agreement  is  better  than  would  be  expected. 

B.  Caliber  . 30 


Data;  w  =  0.012  in 

v  =  2600  ft/sec 
m 

C  =  49.8  grains 

M  =  172  grains 

A  =  0.07355  in2 

D  =  0.306  in 
Chamber  volume  =  0.247  in 


Calculated;  e  =  0.2895 

A  =  0.7973gm/cm3 
k  =  0.94V 
aQ  =  0.0259090 

|  =  0.491 
1/b  =  0.3239 
(RlTe)1/2  =  2211  ft/sec 
B  =  1.157 


*  R.  H.  Kent  and  H.  P,  Hitchcock,  Ballistic  Research  Laboratory.  Report 
No.  171  (original  date  April  5,  1929,  revised  18  January  1940). 


We  obtain 


Ve 

1.157  F  =  2I0O  "  °‘491  (Ve  in  ft/sec)  (43) 

The  following  table  gives  values  of  Ve  from  the  above  report  with 
the  resulting  values  of  F. 


Table  IV 

Gun 

Special  Remarks 

ve 

f t/pec 

F 

Automatic  Rifle 

i 

Gal.  0.30  M  1915 

Middle  Gas  Port 

4iip 

0.94 

Small  Gas  Port 

3710 

0.81 

Rifle  Cal. 0.30 
M1903 

First  Barrel 

390£ 

1 

0.87 

Second  Barrel 

372<j) 

0.81 

Machine  Gun 

Cal. 0.30  M1917 

3839 

1 

0.85 

Mean  = 

0.86 

4*  Discussion  of  the  Results 

Table  V  summarizes  the  mean  values  of  F. 
Table  V 


Gun  Mean  F 


240mm  Howitzer 

Model  of  1918  Ml 

1.00 

3"  Seacoast  Gun 

Model  of  1898 

0.97 

Caliber  0.50 

1.02 

Caliber  0.30 

0.86 

Heglect  of  gas  friction  and  of  the  chambrage  effect  is 
expected  to  be  the  most  serious  in  the  case  of  the  smallest 


calibersf  with  a  corresponding  overestimate  of  the  recoil 
momentum  after  ejection  and  a  resulting  smaller,. vaiueof or . F . 
it  is  therefore  very  reasonable  that  the  value  for  the  caliber 
0.30  comes  out  somewhat  lower  than  the  others.  The  larger  value 
1.02  for  the  caliber  0.50,  however,  is  unexpected,  but  it  is. 
based  on  only  three  rounds  as  compared  with  thirteen  for  the 
caliber  0.30. 


It  is  appropriate  at  this  point  to  say  a  little  more  about 
neglect  of  the  chambrage  effect.  In  small  arms  the  ratio  of  the 
cross-section  of  the  chamber  to  that  of  the  bore  may  amount  to 
as  much  as  2.  It  might  therefore  seem  surprising  that  neglect 
of  the  retarding  force  on  the  joining  cone  does  not  lead  to  very 
serious  error  in  such  cases.  The  answer  is  that  our  assumption 
of  a  uniform  cross-section  equal  to  that  of  the  bore  largely 
compensates  for  this  effect,  since  it  assigns  too  small  a  value 
to  the  area  of  the  breech.  If  there  were  no  pressure  drop  between 
the  breech  and  the  entrance  to  the  bore,  the  compensation  would 
be  exact.  Actually  the  latter  drop  is  only  a  small  fraction  of 
the  total  pressure  drop  from  breech  to  muzile,  so  that  the  com¬ 
pensation,  although  not  complete,  should  be  almost  so. 


For  guns  other  than  small  arms  it  appears  that  the  value 
unity  should  be  satisfactory  for  F.  Thus  we  may  use  Eq.  (40) 
with  no  empirical  correction  factor. 


IV.  fS^immarv 

For  guns  other  than  small  arms  the  effective  velocity  of 
escape  v  of  the  powder  gas  after  ejection; is  given  by  the  following 
equation! 


K  =  Vvm  =  I  +  hf^>  (RlTe)1/2/Vm  <*°> 

V  =  muzzle  velocty 
m 

e  s  ratio  of  charge  weight  to  projectile  weight 

y  =  ratio  of  specific  heats  of  th£  powder  gas 

h  =  1  -  ^  aQ,  where  n  s  l/(y-l )  |  (25.1' 

aQ  is  given  by  calculating  the  following  quantities; 


E  2  2^5+1)  ^ 

*>!  =  -fs  [a  +  f  n  E)1/2  '  (36) 

E*  5  E  -  n(n-l)b|  (37) 


*The  gas  friction  effect  4#  e: 
■Biniy  because  the  small  ea~ 
larger  guns. 


eted  to  be  more  Important  with  small  calibers 
guns  are  longer  In  cal&bers  than  the 


(38) 


* ■  k  E1  +  f n  E’>1/2  - x] 

Then 

ao  =  1+b 


(39) 

(16) 


Rjlg  =  Xi  -  Xf^(1+k^/&)V^  where 


(17) 


k  = 


X  =  specific  force  of  the  powder 

25.8  w/D 

T  1/2 

A^I5%  “  ^ 


(13) 


w  =  web  thickness  of  the  powder 

*  j  l/2 

DjE effective  caliber  *  (“  bore  cross-section) 

^  i  ‘  * 

(w  and  D  are  to  be  expressed  in  the  same  uijiits) 

A  *.  density  of  loading  in  gm/cms 

T  -  adiabatic  flame  temperature  of  the  ponder  in  degrees  Kelvin 
o  ~ 

1  -(±-±) 

Vo 


i/b  = 


2n+3 


Equation  (40)  has  been  shown  to  be  satisfactory  for  large 
caliber  guns  for  values  of  e  as  large  as  1/3.  For  values  of  e 
less  than  0.1  it  is  sufficiently  accurate  to  use  h  =  1  and  0  -  3. 
The  somewhat  lengthy  calcnlfttion  o-f  aQ  can  then  be  omitted.  It 

was  desired  to  obtain  an  equation  which  would  be  Satisfactory  for 
values  of  e  as  large  as  1.  For  such  a  value  of  e,  h  is  about  0.95 

l/2 

and  the  factor (jj^|4f)  does  not  have  to  be  modifiedJ  the 
extension  to  such  Jfalues  of  e  does  not  therefore  involve  more 


than  about  a  %  correction.  Bven  if  the  latter  correction 
itself  should  be  in  error  by  as  much  as  20% ,  the  added  error 
would  be  only  1%.  The  equation  is  therefore  ejected  to  be 
satisfactory  for  values  of  e  as  large  as  unity. 


^John  P.  Vinti 


I 
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APPENDIX:  CALCULATION  OF  aQ  FKQM  £ 


1 


From  (23)  and  (24)  ae  have 

e  -  2(n+l)  a0(l-ao)-n“1j'(l-aof.2),:1d|J 


>,rith  the  definitions 

and 


2(n+l) 

a^ 


£  E 


1-i 


=  b# 


o 


we  have 


N0w 


1 


n 


E  =  b  /  C-^rt — )  dH 


(Al) 

U2) 

(A3) 

(M) 


1  -  a_ua  a^ 

l_a  9  =  1  +  =  1  +  b(l-f;£) 


o 


Thus 


I 


E  =  b/[l  +  b(l-p£)] 
o 


n 


( A5) 


(A6) 


On  expansion  oi  the  integrand  in  powers  of  b  and  integration, 
we  have : 

E=b  [l  +  §nb+  ^n(n-l)b2  +  y^n(n-i)  (n-2)b®+  ^|^n(n-i)  (n-2)  (n-3)tA+„  .  ](£#• 


For  n-4  the  series  terminates  with  the  last  term  written,  and  for 
b  i  6.1,  E  -  Cr.lfl*  +  0*2666  ■+  0,032  +  0*0013,2857  +  0. 0000,40651  = 

0.13005359,  so  that  £  =  10E  =  1.3005359. 

Inspection  shows  that  the  error  involved  in  stopping  with  the  term 
-^r  n(n-l)b3  would  be  only  0.23?.  This  fact  suggests  using  only 
the  first  three  terms  and  solving  the  cubic  for  b. 


cubic 


For  values  of  b  =  0-1,  i.e.  for  &  =  1*3,  vve  may  solve  the 

b  h-  #  hbH  +  -^r  n(n-l)b3  =  E 


(A3) 
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to  high  accuracy  by  successive  approximations. 

In  (18)  neglect  or  the  term  in  b3  gives  as  a  first  approximation 
for  b: 

bi  -  £  [a  ♦  !  n  E)l/2  -l]  U9) 

Insertion  of  b^  into  the  cubic  term  in  (A8)  gives 

b  +  |  n  bz  =  E  -  ^  n(n-l)b®  »|e«  (410) 


Solution  of  (A10)  then  gives:  -  ! 

b  =  -n  [(1+  I  n  E')1/2  -  y 


(ill) 


Ea.  (ill)  completes  the  derivation  of  the  method.  For  n  -  4  an(i 
e  =  1.3005359  the  value  of  b  is  known  to  be  0.10000000.  The 
above  approximate  procedure  gives  b  =  0.09998477  and  a  =0.09089650 
to  be  compared  with  the  accurate  a  -  O.l/l.l  -  0.0909,0909* 


The  error  is  only  1  part  in  7000.  By  Eq.  (33)  the  accurate  arid 
approximate  values  of  aQ  give  respectively  for  I/O  the  values 

0.3010  and  0.3011.  The  above  procedure  is  thus  adequate  for  the 
calculation  of  6  for  values  of  e  -  1. 


> 


It  may  be  asked  why  the  procedure  gives  b  correct  to  0.015# 
when  the  cubic  for  E  gives  an  error  as  large  as  0.2#.  The  answer 
is  clear.  For  a  given  value  of  E,  bi  is  an  overestimate  for  b 
because  of  neglect  of  powers  higher  than  b  .  Use  of  b^  in  the 

cubic  term  thus  overestimates  the  latter,  thereby  compensating 
for  neglect  of  powers  higher  than  the  cubic. 


-21- 


TABLE  OF  SYMBOLS 


SYMBOL  DEFIHIXIOH 

A  Cross-sectional  area  of  the  bore 

a  Dimensionless  parameter  in  Kent’s  special  solution 
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b1  First  approximation  for  b 
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E’  Abbreviation  for  E  -  ^  n(n-l)b^ 
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Empirical  factor  to  be  adjusted  to  give  correct  recoil 
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Mass  of  recoiling  parts 
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SYMBOL  DEFINITION 

I 

p  Breech  pressure  at  ejection 

pe  Space-mean  pressure  at  ejection 
Gas  constant  per  unit  mass 

r -L  Ratio  of  density  at  the  breech  to  spape-mean  density 

,  Ratio  of  pressure  at  the  breech  to  sphce-mean  pressure 
S  The  integral  ^ (l-aopa)n  dp 

„„  T0  Adiabatic  flame  temperature  of  the  ponder  gas  (?K) 

l 

Tg  Space-mean  temperature  of  the  gas  at  1  ejection  (®K) 

t  Time  measured  from  ejection 

Vg  Effective  velocity,  of  escape  of  the  powder  gas 

Vm  Muzzle  velocity  of  the  projectile 

Vpf  Maximum  or  final  velocity  of  free  recoil 

V  Velocity  of  free. recoil  when  the  projictile  base  is  at 

rm  the  muzzle  1 

I 

w  Web  thickness  of  the  powder 

X  Ratio  of  the  total  volume  from  breech  to  muzzle  to  the 
cross-sectional  area  of  the  bore 

Y  Ratio  of  specific  heats  of  the  powder  gas 

A  Density  offloading 

h  l/b  is  the  fraction  of  the  mass  of  the  powder  to  be 
added  to  the  mass  of  the  projectile  fqr  computation 
of  the  total  kinetic  energy  of  projectile  plus, powder 

s  Ratio  of  the  mass  of  the  powder  to  that  of  the  projectile 

X  Specific  force  of  the  powder 

H  Dummy  variable  of  integration  in  the  integral  S 

p  Density  of  the  powder  gas  at  the  breech 

Space-mean  value  of  gas  density  in  the  gun 
pQ  Breech  density  at  ejection 

p_  Space-mean  density  at  ejection  =  C/(A£) 

Time  constant  for  the  rate  of  fall  of  pressure  after 
ejection 

q>  Rate  of  mass  efflux  of  powder  gas  after  ejection 
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